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Abstract 

A mechanism is discussed to obtain light scalar fields from a spontaneously broken continuous sym- 
metry without explicitly breaking it. If there is a continuous manifold of classical vacua in orbit space, 
its tangent directions describe classically massless fields that may acquire mass from perturbations of the 
potential that do not break the symmetry. We consider the simplest possible example, involving a scalar 
field in the adjoint representation of SU(N). We study the scalar mass spectrum and its RG running at 
one-loop level including scalar and pseudoscalar Yukawa couplings to a massive Dirac fermion. 



1 1 Introduction 

In a theory with a global or local symmetry group G spontaneously broken by a multiplet of scalar fields, 
the classical vacuum manifold y is globally invariant under the action of G. The set y/G of classical 
minima in orbit space is often discrete. When y/G is continuous its tangent directions are null vectors 
of the second derivative of the potential and, therefore, describe classically massless modes which arc not 
Goldstone bosons (henceforth referred to as GBs) associated to the spontaneous breaking of G. It may 
happen that these tangent modes are GBs of a larger symmetry group G' D G of the scalar sector, explicitly 
broken by other interaction terms invariant under G but not G'. In other words, they may be pseudo- 
Goldstone bosons (henceforth PGBs) related to G'. If no such larger symmetry G' exists, the tangent modes 
may acquire mass from additional small G-invariant terms in the scalar potential which cause deformations 
of y. We have then light scalar modes without explicitly breaking the symmetries of the theory, neither at 
tree level nor radiatively. 

In this paper we consider the simplest possible model of light PGBs without explicit symmetry breaking, 
involving a scalar field in the adjoint representation of SU(iV), N > 3. We study the mass spectrum of the 
model, and its renormalization when the scalar fields couple to Dirac fermions. In dimensional regularization 
fermion loops may have a large effect on the renormalization group (henceforth RG) running of scalar 
couplings and masses. In particular, as discussed below, fermionic radiative corrections may cause the mass 
of tangent modes to be decreasing functions of the renormalization scale. Because our aim is to investigate 
the mechanism outlined above to obtain light scalar fields, for simplicity we omit gauge interactions. 

In section [2] below we briefly review the extremum analysis for a quartic potential in the adjoint repre- 
sentation of SU(N). In section [3] we consider the classical mass spectrum in the vacua of interest to us. The 
renormalization group evolution of that spectrum is studied in section [5J where we also consider in detail 
the conditions for stability of the vacuum. We try to keep the overlap of these sections with the previous 
literature, particularly [TJ[2], to a minimum. Final remarks are given in section [5j and some related material 
is gathered in two appendices. 
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2 Classical vacua 



We consider a scalar field 4>(x) in the adjoint representation of SU(iV), N > 3. The most general renormal- 
izablc potential Vm ■ su(N) — > M invariant under SU(iV) is, 

Vn(4>) = -^Tr (0 2 ) + cTr (0 3 ) + a(Tr (0 2 )) 2 + 6Tr (0 4 ) . (1) 

We restrict ourselves to the strictly renormalizable case, a =/= =/= b, with spontaneous symmetry breaking, 
m 2 > 0. We introduce the reduced potential Vjv(x) = Vv(diag(xi, . . . , xjv)), with J2f=i x j — 0- Ignoring 

the tracelessness constraint for the moment, the condition for Vjv(x) to be bounded below is that Vjy(x) = 

a(Y^,f = i x 2 ) 2 + bJ2jLi x j > f° r all Since V^\x) is an homogeneous function of x of even degree, this is 

equivalent to requiring that the minimum value of Vjy (x) over the unit sphere be positive. That minimum 
value is a + b if a > and Na + b if a < 0, so a necessary and sufficient condition for Vjv to be bounded 
below is 

a + b > and Na + b > . (2) 

When the constraint Ylj=i x j = is taken into account the conditions (J2J are only sufficient, not necessary. 
We will always assume, nevertheless, that a and b are chosen so that © are satisfied. 

To find the extrema of Vjv (x) we consider the unconstrained extremization of Cn(x, A) = Vjv(xi, . . . , xjv)+ 
X^2^ =1 Xj — 0, with A a Lagrange multiplier. £at(x, A) is a quartic symmetric polynomial in Xj. If x° is 
an extremum oi Cm, N > 4, there cannot be more than three different values among its components if 
6 7^ 0. To see this, we assume otherwise. Since Cm is invariant under permutations of Xj, we can always 
assume that x° is such that x°i < x°2 < x°3 < x°4. Thus, x" is a point interior to the open domain 
{xi < X2 < X3 < X4} C M. N , and within that domain we can use as coordinates s%, s%, S3, S4, X5, . . . , xat, with 
Sfc = Y^ A j=i x j- necessary condition for x° to be an extremum is then {DCm /dsi)\x° = b = 0. Therefore, 
for b 7^ all extrema of the potential must be either of the form, 

x" = ( t?i, . . . , T71 , 772 , t?2 , ??3, ■ - , ^3 ) , ni + n 2 + n 3 — N , ni, 2 ,3 > 1 (3a) 

ni n 2 113 

nirji + n 2 i] 2 + n 3 r/ 3 = , (3b) 



or of the form 



x° = ( 771, . . .,rji , ri2, ■ ■ ■ ,Tfr ) , k 1 +k 2 =N, fe 1)2 > 1 , fciTji + /c 2 ry 2 = . (3c) 

Since Vjv is symmetric in its arguments Xi , any permutation of the components of x° is also an extremum. 
We can then adopt the convention rv\ > n 2 > n 3 . In some exceptional cases discussed below, an extremum 
can be both of type (|5a|) (with, say, 771 =772) and of type (f3"c) . 

In the case of extrema x° of the form ([3a| with 77^ 7^ 77^ if i 7^ 7 = 1, 2, 3), we can argue as above to 
conclude that 77^ must satisfy {8Cm / ds 3 )\ x o = (c + 4/3651)1^0 = or, more explicitly, 

3 c 

m+m + m = -|t ■ ( 4 ) 

Taking into account (f3"a[) , (]3b[) , and (JH) , the extremum equations dC m /dxi = can be reduced to 

TTl 2 9c 2 

(2nm + 6)t/ 2 + (2n 2 a + &)?? 2 + (2n 3 a + 6)77! - — - — = . (5) 

The constrained extrema of Vat(x) are then either of the form (|3a|) . with 771,2,3 determined by (|3bj) . ((4]) 
and ((5]), or of the form (|3c|) . The latter can be obtained by setting 713 = in the solutions of type (l3a|) . 
with eqs. (|3bp and ([5]) determining 771,2 in terms of 773 as defined by (|4|). In order to find the extrema of 
Vjv(x) we consider the case c = first, and the general case afterwards. Since Vjv(x) is invariant under the 
transformation x — > — x, c — > — c, we can restrict ourselves to c > . 
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Even potential (c = 0) The solutions to eqs. (I3b[) . J5]) in this case are, with notation analogous to 
2/°7H,n 2 ,n 3 = (£l>--->£l>£2,---,&,£3,---,£3) , ^i + n 2 + n 3 = N , Ui, 2 ,3 > 1 (6a) 




(6b) 



rai ,712,13 

fim.na.n, = •'•••'• , « + ^ 1 '" 2,n3 6, Q«i.«=.™3 = jv( nina + n2 „3 + „ 3ni ) _ 9nin 2 n 3 , 

nni.na.Ba = ^2 _ + ^ + _ 



Besides ±J/° ni „ 2 n , all permutations of their components are also solutions to the extremum eqs. Clearly, 
if = n 2 = 713 © reduces to the trivial solution. If not all 7ij arc equal it is possible to prove that 
< fl2 un2 ' n3 , < ft™ 1 '™ 2 '" 3 , and 1 < Q£i.n 3 ,n3/nj» 1 ' na ' n3 < ^ provided © holds. Therefore n ni ,„ 2 , n3 > 0, 
and these extrema exist. 

If N = 3n it is possible to have ni = n 2 = n.3 = n. In this case eqs. (|3b[) and ([4]) are not independent so, 
instead of a discrete set, the extrema lie on the curve y° n ,n,n defined by 

&+6+6=o. e 1+ a + e 3 = <, (7) 

If |J2| holds, then (2na + 6) > 0. For iV = 3 these arc minima. 

We can find all two- valued extrema y ni ^ n2 directly by setting n 3 = in y ni ,n 2 = y°m,n 2 ,o- Using © 
it can be shown that there cannot be integers ni.2,3 > 1 and fci )2 > 1 with ni +n 2 + TI3 = N = fci + fc 2 such 
that y°„ lin2; „3 and y°fc li / C2 are equal, up to a permutation of their components. If TV = 3n, however, some 
of the y°n,n,n may have only two different values among their components. In fact, y° 2n n with components 
£i )2 as given by (|6b|) with n$ = are easily seen to satisfy (J7J, with ^ = £1, ^ = £i> £3 = £2- Thus, ±y° 2nn 
and their permutations are part of the set 3^°n,n,n- 

Non-even potential (c > 0) The solutions to the non-homogeneous linear eqs. (|3b[) and (H|) can be 

written as rji = <r£j + pi, i = 1,2,3, with the solution to the homogeneous eqs. given by ([6]), a a free real 
parameter, and pi a particular solution to ([3b| and ((4]). In the case ni > rt 2 > n 3 > 1, pi can be taken to 
be, 

3c n 2 n 3 3c nin 3 3c ri\n 2 . . 

Pi = 77 i P2 = 77 , p 3 = 77 , n y = ni - rij , 8 

4bn 12 n 31 46 ni 2 n 23 40n 23 n 3 i 

and similarly in the cases n 2 = n 3 and n\ = n 2 . There are no solutions with n\ = n 2 = n 3 . The parameter 
a is determined from ([5]), which reduces to 

m 2 9 / , % ^ 1 , , , 9 to 2 9c 2 

— <r 2 + 2J2an* + &)&pj CT + - ^(2a^ + 6)p 2 = — + — . (9) 

»=i j=i 



The extrema of type (|3a[) are then explicitly determined by (|6J), JSJ), and the two solutions to (|9J). There are 
also extrema of the form (|3"c)) wich are given by © and ijHJ), with n 3 = 0, and 



3c Til — 712 

C7 — 2z(7 — 1 = 0, z = / 10 

4to y/il nunai0 



We will denote the two possible values of a by a± — z±\/ z 2 + 1, and the corresponding extrema as xP nitn3 ±. 
Since in this case pi, 2 = 0, we have x° nunat ± = u±y° ni n2 . 

Provided the conditions (J2J) are satisfied, for fixed values of a, b 7^ 0, m 2 > 0, and 774.2,3 _ 1, there exists 
at least one value of c 2 such that two of the rji in ([3a| are equal. For those values of c the extremum x° can 
be written both in the form (f3aj) (with n 3 > 1) and (|3c|) . We will not dwell longer on this issue because it 
will not be of importance in what follows and, furthermore, because such fine tuning of parameters is not 
preserved under renormalization. 
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3 Mass spectrum 



At each extremum x° the Hessian matrix for the potential Vn(x) can be diagonalized (see appendix [AJ to 
obtain the corresponding mass spectrum. In the case c = this diagonalization can be carried out explicitly 
and the extrema y° ni , n2 , n3 and y° „ 2 completely classified. When c^Oa complete classification is more 
difficult to obtain by direct computation. A proof that all extrema x° niin2 na with three different components 
are saddle points is given in [2]. At an extremum x nufl3 ± (n\ > n 2 > 1) the Hessian matrix has the four 
eigenvalues, 

wi = , [1] {2nm 2 + 1} 

w 2 =m 2 (4 + l), [1] {1} 

m 2 N ( + %^o* ± b) , [m-1] {«?-!} (ID 



w 3 



ni - n 2 O ni) „ 2 



uji = m N [ h — a±b , n 2 - 1 {n 2 - 1} . 

\ ni - n 2 Sini,n 2 ,0 / 

The square brackets indicate "reduced" multiplicities, corresponding to the Hessian matrix of the reduced 
potential Vn(x). These are multiplicities in orbit space (see appendix [B| . The "total" multiplicities obtained 
from the Hessian matrix of the full potential Vm{4>) on su(iV), indicated in (fTTT) in curly brackets, can be 
computed from reduced multiplicities f appendix lB) . or simply by noticing that the stability group of x ni ,n 2 ,± 
is SU(ni) x SU(n 2 ) x U(l) if n 2 > 1, and SU(m) x U(l) if n 2 = 1. One null eigenvalue (uj\ in (fTTjl ) is an 
artifact of the projection on the constrained subspace Y2i=i x i = (see appendix [A"|) . The eigenvalues 102,3,4 
give the squared masses of three massive multiplets. The mass-squared average and difference between the 
SU(m) and SU(n 2 ) multiplets is0 

1 , , m 2 N 2 b of Aza+ n-i ~ n 2 \ , , 

^3 + ^) = —- 4, U3-uJi = m 2 N[ —-3^ -&4 • (12) 

Notice that 0J3 + L04 < if b < and n 2 > 1, so x° ni , n2 ,± cannot be minima in that case . In some open 
regions of the plane (a, b) it is possible to choose c = c(a, b) so that 103 — L04 — but, with the exception of 
the case ri\ — n 2 , c = 0, such fine tuning is not preserved under renormalization. 

For even potentials (c = 0) if N = in there is an extremal manifold y° n ,n,n defined by ||7J), containing 
in particular the extrema ±y 2n „• We are interested in the mass spectrum at those extrema and, for c 7^ 0, 
at the related extrema x°2n,n,±- We consider the case N = 3n with n > 1 first. The mass-squared spectrum 
of Vn(x)\ c= o at the extrema satisfying Q and such that & 7^ £j ii i j is (see appendix [A"|) 



2i = , [2]{6n 2 + 2} , 5 2 = 2m 2 , [1]{1} , 



53 = 126(^-^1, 5 4 = 12Me 2 2 --^j , 5 5 = 126k 3 2 - [n-l]{n 2 -l} 



(13) 



where we used the same notation as in (TTT|) and k 2 is defined in 0. All the eigenvalues on the last line 
have the same multiplicity. The null eigenvalue u>\ has two eigenvectors in orbit space, su(A r )/SU(iV), as 
indicated in (fT3"j) . One of them is a spurious projection mode. The other one lies on the subspace tangent 
to the 1-dimensional manifold defined by 0, on which Vn{x)\ c =o is constant. The stability group of these 
extrema is SU(n) 3 x U(l) 2 , so we expect to have 6n 2 GBs plus one spurious massless mode, three massive 
SU(n) multiplets, and two independent massive U(l) modes. Out of the latter, one lies in the tangent 
subspace and remains massless at tree level. Since it is impossible to satisfy (0 and to have (£ 2 — k 2 /6) > 
(or < 0) simultaneously for i = 1, 2, 3, these extrema are saddle points for any values of a and b. 

The extrema of Vn(x)\ c =o satisfying with £1 = £ 2 are ±y° 2nn and their permutations, with stability 
group SU(2n) x SU(n) x U(l). The eigenvalue spectrum at y° 2n ,n can ^> e obtained either from (fTTj) (with 
c = 0, ni = 2n 2 = 2n) or from (fl3|) (with ^1 = = ±k„/V6, ^3 = — 2^i). The result is that there is a null 
eigenvalue with total multiplicity 8n 2 comprising An 2 GBs, one spurious projection mode, and An 2 — 1 PGBs. 
The latter form an SU(2n) multiplet containing the mode lying in the tangent direction to the manifold 



1 The average mass in (1 1 2 I t is slightly different from the result in [2], 
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in orbit space. Since the tangent mode is massless at tree level, the mass of the entire multiplet must vanish. 
The spectrum contains also a U(l) mode of squared mass 2m 2 and an SU(n) multiplet with n 2 — 1 modes of 
mass 66k 2 . We see that for b > all Hessian eigenvalues at ±y 2n n are non-negative. Expanding Vn(x)\ c =o 
in power series about ±y° 2n n through third order, however, shows that these extrema are not minima but 
saddle points. 

When c there are no solutions to ((7]). The only extrema related to the extremal manifold y° n ,n.n are 

x°2n.n,± as defined in with x°2n,n.± ► ^ z V°2n n- As discussed above, we need only consider the case 

c > 0. In this case it is immediate from (jlip that x°2 n .n,— cannot be a minimum, because at it W3 < 0. On 
the other hand, at x°2 n ,n,+ the non-zero eigenvalues in the spectrum (jlip reduce to 

uj 2 = m 2 (l + al ) = 2m 2 + (cm) , u 3 = 6m 2 za+ = -i + O (c 2 ) , 

^ V 2na + (14) 

uj 4 = 6k 2 (2na(l - cr 2 ) + 6) = 6k 2 & + O (cm) , 

with k 2 defined in j7]) and z and cr+ in (|10|) . The counting of modes is the same as in the case of y 2n n 
discussed above, except that in (fT4|) the SU(2n) multiplet has a mass y/uJs = O(^Jcm). The eigenvalue L02 
is obviously positive and, since we are assuming c > and also W3 > 0. The conditions on the coupling 
constants for L04 > 0, and for L04 W3, can be summarized as follows 

=> UJ4 > lu 3 > , 

a>0, b>0 and < A /l:7z=a < £ < \/ 35= 6 "4 > , 




W3 > > UJ4 , 

a> , &<0=^a>3>0>W4, (15) 

' Ana + b < =^ w 4 > w 3 > , 

a<0, 6>0 and ^ 4na + b > and £ < ^-j=M== => uj 4 > oj 3 > , 



4na + 6>0 and «/f .„ b , . < — =4> w 3 > w 4 > . 

We remark that (fT5|) holds under the assumptions AT = 3n, n > 1, c > 0, and ([2]). 

In the case N = 3 (n = 1), the eigenvalues of the Hessian matrix of Vs(a;)| c= o at the extremal manifold 
y° 1,1,1 defined by ||7J) are uj\ — (with reduced and total multiplicities [2]{8}, resp.), and LO2 = 2m 2 (with 
[1]{1}). One of the null eigenvectors in orbit space is a spurious projection mode, and the other one lies on 
the subspace tangent to on which V3(x)| c =o is constant. Since there are no further eigenvalues, these 

extrema are minima independently of the values of a and 6, as long as the potential remains bounded below. 
V3(0)| c= o is invariant under SO(8) (because Tr0 4 = l/2(Tr(/> 2 ) 2 in su(3)), spontaneously broken to SO(7). 
Thus, the 7 massless modes belonging to ui\ are GBs arising from this spontaneous symmetry breaking. At 
the extrema ±y° 2 1 the SU(3) subgroup of stability is SU(2) x U(l), so there are 4 GBs, a massive U(l) 
mode and an SU(2) triplet of PGBs, which can acquire mass from SO(8) breaking interactions. This triplet 
contains the tangent mode. When c ^ the spectrum at x°2,i,+ is as in (|14p . with n = 1, and with uj± 
omitted. The triplet of PGBs acquires a mass of 0(y/cm) from the explicit breaking of SO (8) due to c ^ 0. 

Unlike the case N = 3, for N > 3 the symmetry SO(7V 2 — 1) is explicitly broken by the dimension 4 
operator Tr(0 4 ) in Vn(4>) so we expect it to play a limited role in the theory, and not to have any influence 
on the scalar mass spectrum for \b\ > \a\. 



4 Renormalization group evolution 

We consider now the RG evolution of the masses and the parameters in the potential. We consider the 
interaction of the scalar field <f> with a fermion field in the fundamental representation of SU(-/V), 

L = ^Tr (d^d»<j>) - V N {cj>) +~$iW>- W + iM 575 )V - ^(9 + i9575)i> ■ (16) 
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The one-loop RG equations in MS scheme for the Yukawa couplings and fermion masses are, 
^ = & = CT(^ 2 " 3)55 5 2 + 7rir T (N -1)(N + 3) 5 3 



,f - />* - - ^ + 3 ^ 2 + siv^ 2 - 3 ^ 3 



(17a) 



(17b) 



d/x 8tt 2 7V 
For the dimensionless parameters in V/v we have, 

% = & = + 7)fl2 + ^ (iv2 + 3)52 + (2A?2 ~ 3)afc + i aff2 (18a) 

M~r~ = Pb = 4w^ N2 ~ ^ + + ^ - A(5 2 + flg) 3 , (18b) 
and for the dimensionful ones, 



^ = i € ((" ~ 7?) 1 + 2a + ^" 2 ) ~ 2? (Ms + Ms9s)<! ' 3 +9 »> 1 1 !N " 



- 1™ 2 (<~ 2 + ^ + t) - ^<" 3 " 4)c2 (19b > 



+ £2 (3 2 + 3 5 2 )(^ 2 + M|) + -j 0;M + ff5 M 5 ) 2 



7T^ 



Some comments about these equations are in order. /3 fl and /3 SB have the same sign as g and 35, resp., so \g\ 
and 1 are monotonically increasing functions of (X, as expected. The coupling a is monotonically increasing 
(see below), but \a\ and b need not be and, therefore, the stability of the potential is not guaranteed, as 
discussed in detail below. The sign of the fermion contribution to flb is independent of N, and of whether 
the scalar field is in the adjoint representation of SU(iV) or XJ(N). In the case N = 1 our sign agrees with 
textbook results (see ch. 5 of [3] and ch. 6 of [4]). If M = = M5, then fidc/du oc c, due to chiral symmetry. 
Eqs. Ifl7aj) . [flS]) and lfl9b|) agree with [5] when 35 = = M = M5 = c = On An important verification is 
that for N = 3 the quantities /3 a 4- 1/2/3;,, /idc/dfi and [idm 2 /dfi should not depend on a and b separately, 
but only on a + 6/2. This test is passed by jJSJ and ([15]). || 

The running Yukawa couplings are easily obtained by setting g(t) — gif(t), g${£) — g?af{t), f(0) = 1, 
with gi, g$i the initial values, f(t) an auxiliary function and t = log(/x / /j,q) , with /i /io > and /io a 
reference MS scale. Solving (I17a|) for f(t) we get 

f,A _ 9i M 95i , 4?r 2 iV 

ff(t) -vr^' ^"Tr^W f "(^FWTl' (20) 

The fact that g and g§ enter asymmetrically in /3 S and /3 g 5, and therefore also in if, is due to the scalar wave 
function renormalization not depending on 175 at this order. From (|20|) we see that the ratio g/g§ is RG 
invariant. With (|20|) . (|17bj) can also be integrated 

ft , , »^ ( 1 v 1 , , „ „ , ( 1 x72 



gf + 9ii \i-t/t{J gf + gi \i-t/t {/ 

M 5 (t) = -J^ ( g,U, + g 5i M 5i ) ( -4/r) 1 - ^t^i^iMi - 9i M 5i ) (- — \ 2 , (21) 

3 JV 2 -1 , , , N 

71 = -372 - ^—t { {9i +g 5l ) ■ 



2 Except for the terms of 0(g 4 ), which were considered of higher order in [5]. Here we take into account the full one-loop 
contribution. 

3 But failed by the RG eqs. in [2], 
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We remark that 71 = 3/2 + (D(N~ 2 ) and 72 = — 1/2 + (D(N~ 2 ), independently of the asymptotic behavior of 
Yukawa couplings for large N. A combination with a particularly simple analytical expression is g(t)M(t) + 
gs(t)Ms(t) = (giMi +guMu) (1/(1 — i/if)) 7l+1 / 2 . Clearly, the one-loop approximation involved in (fTTj) only 
holds for \g\, \g$\ <C 1, so the exact solutions (|20p and (|2"Tj) are physically valid only for t <C if. 

We consider now (|18a,|) . For a and g 2 fixed /3 a is a quadratic polynomial in b and /3 a < is possible only 
if b lies between its roots. From the expression for those roots we can show that j3 a < implies that (J2j) does 
not hold. Thus, if Vn is bounded below then (3 a > and a(i) is monotonically increasing. Like g(t) and <7s(f), 
a(i) and b{t) also have a mobile singularity at a finite i, whose location t s > we can choose as a constant 
of integration. (There is a related singularity for t < which will not be of interest to us). For g 2 + g 2 large 
enough the evolution of a and b is dominated by those couplings and t s = if. We assume g 2 ,g 2 ~ 0(H) or 
smaller. Under that assumption we have t s < if and, in fact, over large regions of parameter space i s <C if. 

In the region near the singularity b(t) can be neglected in (118a|) and, with (|2"D|) , it can be integrated to 
give 

ir 2 A tt 2 I 



i(t) 



A 



2 ^ 2 + 7 )(if-i)(l-(^) A ) 2(tf* + 7)t.-t 



(22) 



(N + l)(N-3)g 2 + (N 2 -3)g 2 



5i 



2 



(N-l)(N + 3)g 2 + (N 2 -3)g 5l 



The second equality on the first line holds for i <C if . If g% — we must set if = 00 in (f2"2"| . For 1 6j| < a j the 
expression (j22p gives a good approximation to a(i) also in the physically relevant region i <C t s . Furthermore, 
for N > 3 we see that A ~ 1 and from the initial condition a(0) = dj we obtain Q 

TT 2 

ts ~ 2(A/ 2 + 7H ■ (23) 

If |6j| > |eti| the approximate solution (|22p does not hold. There is a critical value t\, however, such that 
a(i) > |6(i)| for all t\ < i < i s . Thus, for \bi\ > |a.j| we can find an approximate solution that interpolates 
between the regime \b(t)\ > \a(t)\ (0 < t < ii) and the regime \b(t)\ < a(t) where (|2"2"j) holds. 

We are interested in the UV stability of the potential and of the vacuum. By UV stability of the potential 
we mean that ([2]) is satisfied for all < i < t s if it is satisfied at t = 0. Similarly, we say that the vacuum 
is UV stable if it is a minimum of Vn for < i < i s . In the discussion of UV stability that follows the 
relevant regime is \b(t)\ < a(t) so we assume < dj for that purpose (but not in the numerical solutions 
given below). 

Stability of the vacuum requires b > (see (fTSj)'). However, the condition 6, > is not enough to 
guarantee that bit) > for all < t < t s , since the last term in (118b|) can drive b(t) to negative values. 
Once b becomes negative it diverges to — 00 at i = t s driven by the last three terms in (|18bp . There is a 
minimal value 6 m ; n > such that b(t) > for all < i < t s if and only if b^ > b min . In order to obtain an 
approximate expression for & m ; n the relevant regime is b < a, so we can neglect the first term on the rhs of 
(|18b[) and integrate the resulting linear eq. with a(t) from (f2"2"| and g, g$ from (j2"D|) , to get, 

This approximate solution, together with (123|) . leads to 



1 (?l±?lif (, _- g f 7(2(^ + 7) 



9l 



Ql) ) . (25) 



The condition &j > 6 m ; n is not very restrictive since, for g 2 ,g 2 i ~ 0(ai), typically b min -c aj. The accuracy 
of the approximate relation (|25p is best illustrated by quoting some representative numerical values. For 
a,i = 1/600, gi = 1/10 = <75j, ([25| gives a result for 6 m i n that differs from the result obtained numerically 



4 From and t s < t f <^=> (iV 2 - 3)^ + (AT - 1)(N + 3)g 2 < Sn 2 N(N 2 + 7)a, . Notice that the rhs is cubic in N and 

the lhs is quadratic. With Sn 2 ~ 80 and g 2 , g 2 ^ ~ 0(ai), the inequality is satisfied in large open regions of parameter space. 
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from (|18b[) by 35% for N = 3, 10% for A = 6 and 3% for A — 12. The accuracy is higher for larger values 
of N. 

Since we require b > for < t < t s , and since a is monotonically increasing, if a$ > then {5} are 
trivially satished for all < t < t s and Vn is UV stable. If a, < 0, but such that © are satisfied at t = 
with 6^ > 6 m ; n , (J2j) can in principle be violated at t > if 6(t) > decreases rapidly enough. This imposes 
bounds on the possible values of a, < 0, gi and g^i, which must be satisfied in order to ensure that A a + b 
remains positive throughout the evolution. Since, however, a(t) is monotonically increasing and b(t) > 0, 
([2|) will be satisfied for all values of t (< t s ) larger than a certain critical value. We will not dwell longer on 
this, and from now on we assume a.; > 0, bi > 6 m in > so that a(t), b(t) > for all < t < t s . 

We assume Vn to be of the form (fT]) with c > 0, the case c < can be obtained by means of the 
transformation <f) — > —<f). The fermion masses M, M 5 and Yukawa couplings <?, g§ can be ^ 0, their signs 
entering the scalar sector RG eqs. (TTS)) . (flT))) only through Mg + M 5 g 5 in (|19ap . With the assumptions 
a.i > 0, bi > b min , and Cj > 0, if Mg + M 5 g 5 < then dc/dt > and c(t) is monotonically increasing, 
therefore positive throughout the evolution. If Mg + M§g$ > 0, however, even if a > it may happen that 
c(t) becomes negative, thus diverging to — oo at t = t s . The classical extremum x°2n.n.+ is a saddle-point 
for c < 0. As in the case of b(t) discussed above, we will have c(t) > for all < t < t s only if c, > c m i n for 
some critical value c m i n . An approximate solution to (|19a|) . valid for bi < aj, can be obtained by substituting 
([22]). p4|) (with exp(g?t/ir 2 ) ~ 1), (JUJ) and (21]) into (fl9a|) to obtain, 

x x 6/(7V 2 + 7) 



\i t/i s / ^6) 
3 2 6 N 2 -6A 2 + 7, _ l 
^2 ft + ^ jv jv 2 -5 ' ' Cl = 27^ ( hm + • U ""-' y "" :; - / - + 554) ■ 



From here, we get, 



Cmin ~ (1 - e-"* B ) Ci . (27) 

Notice that for Mg + Msgs < we have c, > > c m i n , so c(i) > for all < t < t s . Substituting (T2^|) in 
(|27[) we obtain an approximate expression for c m i n solely in terms of initial values. To give an estimate of 
the accuracy of that approximation, for a,i = 1/300, bi — 1/450, gi — —1/15 = g$i, Mi — M^i = rrii — 1, the 
value of c m i n obtained with (|27j) and (|23j) differs from the more accurate value computed numerically from 
OH} by 35%, 23% and 11% for N = 3,6 and 12, resp. 

The evolution of m 2 {t) may be dominated by the term oc c 2 in (|19bj) . for c 2 (t) sufficiently large, causing 
m 2 (t) to decrease and eventually become negative^ An upper bound c max on the initial value Ci exists, so 
that Ci < c max ensures that m 2 (t) > for < t < t s . Using the exact solutions for the fermion parameters 
and the approximate ones for a, b, c given above, we find the following approximate form from m 2 (t), valid 
for < bi < a,i, 



1 \TO / / lt\\ 2 2A 2 -3, gf 
m 2 (t)~[- — mfe^ + mfq 1- -- , r, = — — b t 1 ; ' 



i-t/tj \ 1 v 2< s yy ' ' 7T 2 a 2tt 



= — ( " 9 ^^ c ? + + ft 2 ,)( M2 + M !i) + 2 (ft^ + 5 5l M 5 ^ 2 



1 / A 2 -4 

From this approximate solution we obtain, using ([23 
N 



(28) 




(A 2 - 4) 



(4(A 2 + 7)a t m 2 + (g 2 + g 2 5l )(M 2 + M 2 ) + 2{g l M l + g 5t M 5l ) 2 ) 1/2 . (29) 



For the set of parameters quoted below ([27]) the values for c max obtained from (|29|) differ from those obtained 
numerically from (fT9b|) by 23%, 9% and 3% for A = 3, 6 and 12, resp. 



5 Typically, 0^4 changes sign at lower scales than m? does (see below), in principle precluding the possibility of radiative 
symmetry breaking, at least in this vacuum. 
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(a) (b) (c) 
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Figure 1: RG evolution of m 2 (t)/m 2 and u^2,3,i/Tn 2 , for Mg + M 5 g 5 < 0. In all cases N = 12, aj = 1/300. 
(a) h = 1/300, c t = 1/1000, 9l = -1/20 = M t = 1 = M 5i , (b) h = 1/450, c, = 1/100, 9l = -1/10 = g 5i , 
Mi = 20, M 5i = 35, (c) h = 1/900, c, = 1/200, g l = -1/10 = g 5l , Mi = 3/2 = M 5i . The end values of a 
are: (a) a(8) = 0.09, (b) a(8.5) = 0.11, (c) a(9) = 0.13. 

A further condition for the stability of the vacuum is 0J4 > 0. From (fT5|) we see that such condition can 
be formulated as r(t) < 1, with r — \J N/2a 1 / 2 c/ (bm). A necessary condition for r(t) < 1 throughout the 
evolution is obviously r(0) < 1, 

1 /2 

a/ a < b t rrii . (30) 

The dependence of r(t) on t is not monotonic in general, so imposing ([3D]) and r(t s ~) < 1 is not sufficient 
to guarantee that r(t) < 1 for all t < t s . Furthermore, r(t) may have rapid variations, which makes it 
difficult to obtain approximate expressions for it that remain accurate over broad regions of initial values 
in parameter space. For these reasons, we will not attempt to give a single bound on initial conditions that 
ensures maxo<t<t s {r(t)) < 1. Rather, we require initial conditions to obey (|30p . and study the evolution 
of LV4 numerically. The stability condition 0)4 > implies an upper bound on which is typically more 
restrictive than ([2"9"| . 

Fig. Q] shows the evolution of the mass-squared eigenvalues (Ti"4"|) for Mg + M555 < 0. As shown there, the 
location t s of the RG flow singularity has a small dependence on 6^, not taken into account in (|23[) which is 
strictly valid only for 6j <C a,*. remain much smaller than 0J2 and m 2 throughout the evolution, even for 
b > a, when the broken SO(N 2 — 1) symmetry should be irrelevant. 

Depending on the values of bi and Cj, W3 can be smaller, larger, or approximately equal to uj^ throughout 
the evolution, as seen in fig. [1] and (fTS")) . The effect of fermion masses on the evolution of scalar ones is 
apparent in fig. [TJd. In that figure the value of c is close to the upper bound imposed by (|30p . which is 
reflected in the small value of 074 at the beginning of the evolution. Similarly, a slightly larger value of c in 
fig. would cause the dip at the end of the evolution of W4 to reach negative values. 

For Mg + M^g^ > the RG evolution of the mass-squared spectrum remains qualitatively the same 
as in fig. [TJ except for the case c ~ c min in which ujj,(t) (as well as c{t)) is monotonically decreasing over 
essentially all of the interval < t < t s . This decrease is shown in fig. [21 where we chose a small value for a, 
and bi ~ lOa^ in order to obtain a larger t s and a small value of c m i n . For the parameters used in the figure, 
with C, ~ c m i n , 0J3 is separated from the rest of the spectrum by a factor 10 4 -10 6 throughout the evolution. 

5 Final remarks 

At the beginning of sect. [5] we describe the qualitative motivations for this paper. The two possibilities 
mentioned there are realized in the model defined by |T]) and pH[) . at the classical vacuum x°2n,n,+ defined 
in sect. [3] First, in the case N = 3 the scalar sector is invariant under SO(8), spontaneously broken to SO(7) 
and explicitly broken by the fermion sector and also, softly, by the cubic scalar self-coupling if c =/= 0. That 
explicit breaking gives mass to a triplet of PGBs, as described in sect. [4] When c = there is a manifold of 
minima in orbit space, y°i,i,i described by ([7]), which collapses to the two points x°2,i,± for c^0. Second, 
in the case N = 3n > 3 there is no larger symmetry group containing SU(A r ), and no explicit symmetry 
breaking. When c = the manifold y° n ,n,n of 10 in orbit space comprises only saddle points, and collapses 
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Figure 2: As in fig. [TJ but for Mg + M 5 g 5 > 0. N = 12, a, = 1/1000, 6, = 1/100, c, ; = 5.35 x 10~ 6 , 
gi = 1/30 = g u , Mi = 1/10 = Mm- The end value of a is o(10.5) = 0.05. 

into 

2-°2n,n,± for c > 0. At the minimum ^°2n.n,+ 

there is an SU(2n) multiplet, containing the mode tangent 
to y° n ,n,n when c — > 0, with mass oc y/cm. This happens even if (at a given renormalization scale //) we set 
6 ^> |a| > so that, unlike the case N = 3, there is no hint of SO(A 2 — 1) symmetry in C We remark that, 
as long as b ^= 0, there is no eigenvalue proportional to c in the spectrum (jlip at any extremum x niinjl ± 
except for ni = 2tt,2. 

In section@]we discuss in detail the RG running of couplings and masses. Approximate bounds are found 
on the initial values of b and c which are necessary for UV stability of the vacuum. Numerical study of the 
evolution of to 2 and 0^2,3,4 (figs. [T] and shows that W3 (and, for b <C a, also L04) remains much smaller than 
to 2 . For gM + (75M5 < 0, c(t) is monotonically increasing but for gM + 35M5 > we have c > (as needed 
for U3 > 0) only if c(0) > c m ; n , with c mul approximately given (|27|) . Interestingly, for c(0) close enough to 
(but still larger than) its lower bound the mass y/uJ^ of the "tangent mode" multiplet is a monotonically 
decreasing function of the renormalization scale throughout the evolution. 

For TV > 3 the vacuum x°2n,n,+ is not a global minimum, which in the simple model discussed here is not 
important. In fact, even in realistic models a metastable vacuum is not in conflict with phenomenology if 
long-lived enough. Although the mass splittings we obtain with this model are somewhat modest, we may 
speculate that the mechanism described in this paper could lead to more elaborate theories in which the 
massless modes in the spectrum are removed by the Higgs mechanism and the heavier ones are pushed up 
by radiative corrections to very high mass scales, so that only the lighter massive modes (corresponding to 
LU3 in fig. [2]) remain in the physically accesible spectrum. 
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A Projected Hessian matrix for constrained extrema 

Let V : l w -» R and g a : R N -> E, 1 < a < A < N, be smooth functions. Let Q C R N be the implicit 
manifold defined by the set of equations g a {x) — 0. We assume that rank(9j<7 a (a;)) = A for all x € Q, 
so that Q is a smooth manifold. We are interested in the extrema of V\g : Q — ► R. Let x° be such an 
extremum, c : [—1,1] — > Q a smooth curve such that c(0) = x", and V(t) = V(c(t)). We must have 
V"(0) = diV(x°)ci(Q) = 0, and since this must be true for any vector c(0) tangent to Q at x°, we conclude 
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that diV(x°) — y^t—i ^° a di9a(x°) for some set of numbers A°. Thus, we are led to an extremization problem 
for C(x, X a ) = V(x) — X)a=i ^adig a (x), which is the method of Lagrange multipliers. 

The sign of V(0) — df 3 V(x°)ci(0)cj(0) + diV (x°)ci(0) determines the nature of the extremum x°. Twice 
differentiating the relation g(c(t)) — and expressing diV(x°) in terms of Lagrange multipliers we get 

V(0) = ^V(x°) - E A °« 9 «^°)) ^(0)^(0) . (31) 

For a; to be a minimum (resp. maximum) the expression on the l.h.s. must be positive (negative) for all 
tangent vectors c(0). Therefore, the Hessian matrix in which we are interested is 

H v (x°) = Pg(x°) (h v (x°) - £ KHgA*")^ Po(*°) > (32) 

where (Hv(x°))ij = dfjV(x°) and similarly H ga {x°), and Pg{x°) is the projector onto the subspace tangent 
to Q at a; . In the simplest case in which the constraints g a are mutually orthogonal (dkg a dk9b oc 5 a b) we 
have 

= h - g ^4^) ^ a(x ° )9igg(x0) ■ (33) 

Due to the projection Pg(x°) the matrix Hv(x°) has A spurious null eigenvectors that span the subspace 
normal to Q at x° and can be taken to be dg a (x°). The remaining N — A genuine eigenvectors must be 
orthogonal to dg a (x°), and therefore span the tangent subspace at x°. Their eigenvalues define the nature 
of a; as an extremum of V\g. 

The extremization problems considered in this paper involve a single linear constraint with dg(x°) = 
(1,1,.. .,1). At a two-component extremum x° of the form (|3"c)) , such constraint leads to Hessian matrices 
of the form, 



H V {X°) = 







«; 3 G("2,nij 


U)J n2 + K 4 G^' n ^ 



(34) 



with (Jl™ 1 '™ 2 ) € R™ lX " 2 , (G^ ni,n2 ^)ij = 1, I n the n x n identity matrix, and ^3,4, K1...4 some numerical 
coefficients. Such matrices are easily diagonalized. Let {v[ n \ . . . , v^} be an orthogonal basis for K™ such 
that = (1,1,..., 1). We define the n x n orthogonal matrix 



uffi = — ={v[ n ^)i (no summation over 7) . (35) 



(n) (n) 



The matrix G^ 1 ' 71 ^ = v[ ni) <g> v[ n2) then satisfies (u^ f G^^U^Y = v^nj^i^i- Let U be the 
N x N orthogonal matrix 

u = { u 7\ul>) • < 36 > 

With this definition the matrix Hv{x°)U is block diagonal, with one block equal to W3/ ni _i, another 
equal to u)Al n ,^i, and a 2 x 2 block ( " 3+ Kini K2 Y' 11 ' 12 ] with eigenvalues wi o, one of which should vanish 

due to the projection. At an extremum of type (l3aj) the form of Hv{x°) is analogous to (j34|) . with 9 blocks 
instead of 4, and the diagonalization procedure is the same as above. 



B Mass spectrum in su(7V) 

In the main text and in appendix [Al we considered the extrema of Vn(x) and its Hessian matrix. From 
those we can immediately obtain the extrema and mass spectrum in orbit space. (The latter being the 
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orbifold su(iV)/SU(JV) = Ef /V N , with = {x e R N / J^f=i x j = °l and the permutation group of N 
elements.) In this appendix we extend the Hessian matrix of Vjv(x) to Vn{4>) over su(iV). In order to do so 
we expand 

AT TV N 

4> = Y,^ Eiii)+ E P^^=(e^ + E^)+ 1 Y, ^-^(EM-EW) (37) 

i=l i<j=l * ' i<j=l * 

with (S (t5 ^) mn = S im 5 jn , 1 < i,j,m,n < N, and use {^i}^ U {pij}^ j= i U {%-}^- = i as our coordinates, 
constrained by X^iLi ^ = 00 The Hessian matrix of V/v(</>) at a diagonal matrix ^ = diag(xi, . . . ,xn) is 
easily computed if we disregard the tracelessness constraint, 

8 2 Vn 

——^-(tp) = (H v (x l7 ...,x N )) lk , (38a) 
ocpiocpk 

J* V f (f) = -f V f if) ^ S ikSji , (38b) 

d 2 V N 8 2 V N d 2 V N 

i ( P) = -z — kj-(<p) = ^, — 5 — 0) = ° ( 38c ) 



dpijd(j)k dritjd&k dp i3 dr]ki 

where Hy (x) is the Hessian matrix of Vat (x) ■ The effect of taking the constraint into account is to substitute 
Hv{x) for Hv(x) in p8|) , as discussed in appendix [A] 

The matrix (|38[) is block diagonal. At an extremum of type (|3c|) the eigenvalues of the block Hv(x) are 
^1,2,3,4 as given in (fTTj) . The other two diagonal blocks are already diagonal, as indicated in (|38b[) . Explicit 
calculation shows that the diagonal entries in those blocks vanish if the index (k,l) is such that x"k ^ x°i 
(see (pel) ). If (k,l) is such that x°k = x"i = f?i (resp. 772) then d 2 V / d pud pu = ^3 (resp. L04). Counting 
the number of such pairs of indices gives the multiplicity of in each of the blocks d 2 V/dpkidpki and 
d 2 V I 'drjkidrjki- Adding those multiplicities to the reduced ones in (TTTj) results in the total multiplicities 
shown in that equation. 



6 Alternatively we could use IV— 1 unconstrained coordinates for the diagonal elements, by expanding in a basis with diagonal 
matrices satisfying TrH^ = 0, instead of E^"\ e.g., a Cartan-Weyl basis. 
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